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1.  INTRODUCTION 


We  report  here  on  the  progress  made  in  our  work  on  the  thermal 
properties  of  fluids  whose  molecules  interact  via  two-tody 
non-spherical  potentials. 

To  avoid  unnecessary  repetition  we  refer  the  reader  to  our  1983 
annual  report^ ^  and  concentrate  here  on  the  new  developments  that 
took  place  this  current  year. 

During  the  present  period  we  directed  our  efforts  to  extending 

F21 

the  spherical ization  procedure  of  the  median  potential  to 
calculate  the  thermodynamics  of  fluid  mixtures  of  rigid  diatomic 
molecules. 

Our  approach  in  considering  this  type  of  fluid  mixtures  was  to 

obtain  an  equivalent  single  component  spherical  potential.  In  order 

to  do  that  we  explored  and  generalized  the  van  der  Waals  one  fluid 

[3] 

{ vdWlf )  mixing  rules,  originally  introduced  to  treat  atomic 
mixtures  interacting  via  conformal  potentials.  We  describe  this 
subject  in  Section  2.  The  mixing  rules  were  used  in  conjunction  with 
the  modification  of  the  median  procedure  that  MacGowan  developed  for 

r  4  -] 

Lennard-Jones  (LJ)  homonuclear  diatomic  fluids,  which  he  calls 
r  51 

MED  (12,6).  In  Section  3  we  give  an  account  of  this  "mixing 
rules  plus  medianization"  procedure  as  applied  to  mixtures  of  diatomic 
homonuclear  LJ  molecules,  for  which  computer  simulations  were 
available  for  comparison.  In  Section  4  we  offer  our  conclusions. 


2.  VAN  DER  WAALS  ONE-FLUID  MIXING  RULES  AND  GENERALIZATIONS 


The  vdWlf  model  is  a  procedure  that  uses  a  simplified  potential 
to  represent  the  mixture  interatomic  potentials.  It  was  derived  to  be 
applied  to  the  case  of  conformal  interatomic  potentials,  i.e., 


t 


ij 


ij 


<t>  (r/b  . .) 

*  J 


(1) 


IP 

where  $  is  a  universal  function  (for  example  4[£  -  S  ]  for  LJ 

atoms)  and  e..  and  b..  define  the  well  depth  and  the  separation  at 

1  J  ■  J 

the  minimum  of  4>..  respectively  for  the  interaction  between  species 

*  J 

i  and  j. 


The  prescription  of  the  vdWlf  approximation  consists  of  defining 
an  equivalent  single  component  fluid  of  potential 

rfx(r)  =  ex  4(r/bx)  (2) 


where  the  x  explicitly  indicates  the  dependence  on  mixture  composition, 
and  ex  and  bx  are  given  by 


bx  ■  2>I  -J  b?j 

sbM>,  xj  su 


(3) 

(4) 


where  x^  are  the  fractional  concentration  of  the  components.  This 

equivalent  potential  is  then  used  to  calculate  the  thermodynamic  of 

the  mixture.  The  results  one  obtains  are  very  good  if  the  spread  of 

values  of  the  e.  .'s  and  b..'s  is  not  too  large.  Indeed  the  vdWlf 
ij  [6] 

theory  does  much  better  than  much  more  elaborate  techniques. 


In  the  case  of  molecular  mixtures  of  non-spherical  rigid 
molecules  our  goal  is  to  jointly  utilize  mixing  rules  of  the  sort  of 
the  vdWlf  theory  and  the  spherical ization  procedure  given  by  the 
median  and  its  extension.  The  objective  is  then  to  obtain  a  single 
spherical  "medianized"  concentration  dependent  potential  representing, 
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for  the  purpose  of  calculating  the  thermodynamics,  the  mixture.  This 
has  certainly  been  our  motivation  in  considering  the  mixing  rules  that 
are  the  object  of  this  description. 

For  the  particular  case  of  binary  mixtures  of  diatomic 

homonuclear  LJ  molecules  described  in  deference  1,  one  needs  eight 

parameters  to  describe  the  1-1,  1-2  and  2-2  interactions;  they  are 

£..  =  e..,  b..  =  b...  and  1.,  i  =  1,2,  which  define  entirely 
ij  ji  ij  ji  i 

the  LJ  site-site  (atom-atom)  potentials  of  the  two  species. 

One  has  at  least  two  ways  to  proceed,  either:  ( i )  "to 
medianize"  the  1-1,  1-2  and  2-2  potentials  using  the  MED  (12,6)  and 
then  to  construct  mixing  rules  for  the  three  resulting  spherical ized 
potentials;  or  (ii)  to  define  mixing  rules  to  obtain  ex,  b^  and 
1^  which  define  a  single  component  diatomic  homonuclear  LJ  molecular 
fluid  and  the  "medianize"  it  using  MED  (12,6). 

Clearly,  the  results  one  would  obtain  by  using  (i)  or  (ii)  will 
not  be  the  same.  The  advantage  of  procedure  (i)  is  that,  at  least  in 
principle,  one  could  extend  this  type  of  reasoning  to  treat  mixtures 
of  molecules  other  than  homonuclear  diatomics.  That  one  can  do 
provided  one  has  a  well  defined  medianization  technique  for  such 
potentials.  One  should  bear  in  mind  that  this  is  not  an  obvious  task 
when  the  molecules  do  not  possess  a  center  of  symmetry,  and  it  has 
only  been  done  for  diatomic  non-homonuclear  hard  dumbells.^ 


Even  in  the  simple  case  of  homonuclear  diatomics  LJ  mixture  we 
are  describing,  in  which  the  site-site  potentials  are  conformal,  the 
resulting  medianized  ones  are  not  of  the  same  shape.  Therefore,  one 
must  deal  in  general  with  non-conformal  spherical  potentials,  i.e., 

*ij  =  eij  *(r/bij’  {5) 

where  e..  and  b. .  have  the  same  meaning  as  before  as  a.,  is  a  set  of 
ij  ij  3  ij 

parameters  introduced  to  characterize  the  difference  in  shapes  of  the 
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An  important  realization  of  Eq.  (5)  is  the  exponential  -6 
(exp-6)  potential  where  there  is  only  a  single  parameter  <*..  for 

*  J 

each  <S . ,  and 

'  J 

<*>(£, a)  =  6  exp[a(l-t;)  -  a^‘6]/(a-6). 


[8] 

Ree  J  has  recently  found  that  for  mixtures  of  exp-6  atoms  -  he 
modeled  an  Hg-H^  mixture  at  high  temperatures  and  densities  -,  the 
mixing  rules  of  Eqs.  (3)  and  (4)  and  the  following  extra  mixing  rule 
for  ax  is  very  successful: 


h3 

VxbK 


■  £ x,-  xj 


With  MacGowan  and  Lebowitz  we  researched  the  possibility  of 
making  the  formulation  of  mixing  rules  more  systematic.  We  found  that 
for  the  case  of  conformal  atomic  LJ  potentials  we  were  able  to  clarify 
how  the  rules  of  Eqs.  (3)  and  (4)  come  about  by  use  of  approximate 
scaling  arguments  on  the  radial  distribution  functions  g..  in  the 

'  J 

compressibility,  energy  and  virial  equations.  For  the  case  of  exp-6 
potentials  by  expanding  to  second  order  the  potential  around  its 
minimum  in  a  Taylor  series  and  employing  again  approximate  scaling 
arguments  on  the  g^.’s,  we  show  how  one  can  generate  various 
alternative  mixing  rules  for  ex  bx  and  ax-  One  such  set 
produces  thermodynamic  results  of  accuracy  equivalent  to  that  of  Ree's 
work.  The  details  of  this  work  are  contained  in  Appendix  A,  which  is 
a  preprint  of  an  article  submitted  and  accepted  for  publication  in 
Chemical  Physics  Letters. 


This  is  the  present  status  on  this  approach.  Up  to  now  to  have 
not  applied  this  procedure  ( i )  to  molecular  mixtures.  Before  leaving 
this  section  we  remark  that  our  work  on  mixing  rules  for  non-conformal 
spherical  potentials  is  a  contribution  in  its  own  right  for  such  fluid 
mixtures  (as  in  the  case  of  the  exp-65,  and  not  just  an  intermediate 
step  in  "procedure  ( i ) "  for  molecular  mixtures.  We  explain  procedure 
(ii)  next. 
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3.  A  ONE  MOLECULAR  FLUID  APPROXIMATION  FOR  DIATOMIC  FLUID  MIXTURES 


We  considered  a  mixture,  treated  as  rigid 

homonuclear  diatomic  LJ  molecules,  for  which  computer  simulations  were 

m 

available.  Because  it  happens  in  this  case  that  the  bond  lengths 
are  very  similar,  1 ^  ~  1 2  ~  2.36  A,  it  occurred  to  us  to  apply 
procedure  (ii)  as  defined  in  Section  2.  That  is:  We  first  find  ex> 
bx  abd  1  which  define  a  single  component  equivalent  homonuclear 
diatomic  LJ  fluid  and  then  apply  to  it  MED  (12,6).  For  ex  and  bx 
we  used  the  mixing  rules  of  Eqs.  (3)  and  (4).  For  lx  any  reasonable 
rule  would  give  essentially  identical  results,  we  chose 

Vl/i'i  l7) 

The  MED  (12,6)  d  (r)  potential  was  utilized  then  in  Ross's 
x  CIO] 

variational  procedure  for  obtaining  the  free  energy  and  from  it 
the  pressure  and  internal  energies  of  the  mixture.  The  results 
compare  favorably  with  those  of  the  computer  simulations,  which  is 
remarkable  considering  the  calculational  simplicity  of  our 
formulation.  (It  takes  1  second  of  CPU  time  of  a  DEC20  computer  per 
thermodynamic  p,T  point.) 


We  also  applied  this  technique  to  equimolar  binary  mixtures  of 
rigid  homonuclear  diatomic  mixtures,  of  molecules  denoted  as  N.N,  N-N 
and  N— N  by  the  researchers  that  did  molecular  dynamics  calculations 
employing  them.^1^  They  all  have  e  • .  =  e.,  ,  b . .  =  bw  , 

IJ  Np  1J  Np 

i , j  =  1,2,  where  and  b^  are  the  parameters  used  to 

represent  N2.  In  tflis  nomenclature  N-N  is  identical  to  N2  and  has 

an  elongation  l/o  of  0.3292  while  N.N  and  N--N  have  half  and  twice 

that  elongation,  respectively.  Here  c  is  the  point  at  which  the  LJ 

1/6 

potential  crosses  the  axis  and  is  related  to  b  by  b  =  2  c.  The 


mixing  rules  for  ex  and  bx  obviously  give  e^  and  b^, 

while  for  this  equimolar  mixtures  Eq.  (6)  gi^es  the  irithmetic  mean  of 

lj  and  I2.  The  sparcity,  and  some  errors  we  believe  to  exist  in 

the  reported  data  by  the  authors  of  Reference  11,  do  not  permit  us  to 

draw  definite  conclusions  in  this  case. 


We  also  developed  an  alternative  mixing  rule  to  that  of  Eq.  (6) 
using  the  systematic  approach  described  in  Appendix  A  for  atomic 
systems,  and  in  Appendix  B  for  the  system  under  consideration.  By 
employing  arguments  of  approximate  scaling  for  the  g . . ' s  in  the 

'  J 

virial  equation  for  these  diatomic  mixtures  we  obtained 


e  b^x  =  x .  x .  £  . .  b? .  x . . 

xxx  i  j  ij  ij  ij 


where  x..  e  (1  +l.)/2b... 

I J  •  J  *  J 


(8) 


In  the  cases  reported  in  this  section  the  values  of  1 
obtained  with  Eq.  (8)  are  either  identical  or  almost  identical  to  that 
of  Eq.  (6).  To  text  the  mixing  rule  of  Eq.  (8)  one  needs  to  compare 
its  results  with  computer  simulations  with  1^  i  1^  and  x^  i  x^ 
which  are  not  presently  available. 


We  include  in  Appendix  B  a  detailed  account  of  the  work 
described  in  this  section,  which  is  a  preprint  of  a  paper  with 
Lebowitz  and  MacGowan  scheduled  for  publication  in  the  December  15, 
1984  issue  of  Journal  of  Chemical  Physics. 
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4.  CONCLUSIONS 


Very  significant  progress  has  been  made  in  extending  the  idea  of 
the  median  potential  to  obtain  the  fluid  thermodynamic  property 
mixtures  of  rigid  diatomic  homonculear  molecules.  Also  we  have  made  a 
contribution  in  clarifying  and  generalizing  mixing  rules  for  atomic 
and  molecular  mixtures. 

We  expect  to  be  able  to  continue  to  make  significant  advances  in 
treating  mixtures  of  diatomic  and  monoatomic  species  as  well  as 
mixtures  containing  other  simply  shaped  rigid  symmetrical  molecules, 
such  as  those  of  tetrahedral  structure. 
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VAN  DER  WALLS  ONE-FLUID  THEORY:  JUSTIFICATION  AND  GENERALIZATION 
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Abstract 

We  investigate  a  one  corriponent  molecular  fluic  approx  irr.  at  icn  tor 
conformally  similar  molecules.  Ue  test  this  scheme  cn  (two)  mixtures 
cf  rigid  hcmcnuclear  diatomic  Lennarc-Jones  (LJ)  fluids  for  which  a 
limited  amount  of  information  from  molecular  cynamics  simulations  is 
available.  For  two  components  cf  approximately  equal  bone  length  but 
cifferent  LJ  parameters  cur  results  compare  favorably  with  the  machine 
computations.  From  the  very  few  simulation  data  available  fer 
eauimolar  mixtures  of  molecules  differing  only  in  their  bone  lengths 
we  cannot  reach  any  firm  conclusion.  Alternative  procedures  for 
treating  general  molecular  fluid  mixtures  are  discussec. 

1 .  I ntrccuct icn 

The  development  cf  simple  accurate  approximation  schemes  for  tre 
thermodynamics  anc  structure  cf  dense  fluids  is  a  problem  cf  great 
practical  in  pcrtar-ce  anc  also  a  theoretical  challenge.  While  tre 
e*  isterce  cf  such  appro  irraticr.s  for  fluics  ir  general  1  .  rc  means 
i  prior'  ctvicus  the  results  cn  spherical  single  ccnponent  fluids  arc 
"core"  fixtures  1’  strcrgly  suggest  the  possibility.  ’"his  motivates 


*-  '•*  r-nr  -caress :  Research  School  of  Chemistry,  Australian  ‘.at  icra  1 
r  •  to  ,  IT'l  lcx  A  ,  -  ar  terra  Al  T  Z  6  (.  1 ,  Australia. 
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A  ONE  MOLECULAR  FLUID  APPROXIMATION  FOR  DIATOMIC  FLUID  MIXTURES 
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Table 


fli xture 

Equ i va ) ent 

pure  f 1 u i 

ds 

T  (K) 

V  (cm3 /mo  1 ) 

X 

i 

Ree 

l 

1  1 

1  1  1 

50 

20.0 

0.50 

0.047 

-0.755 

0.048 

-0.736 

.0.048 

-0.737 

0.049 

-0.725 

0.055 

-0.669 

100 

14.0 

0.50 

0.338 

-0-35*4 

0-337 

-0.33** 

0.336 

-0.336 

0.339 

-0.318 

0.365 

-0.202 

O 

O 

fO 

10.0 

0.50 

1.856 

2.79 

1 .858 
2.91 

1 .850 
2.90 

1  .872 
2-9*4 

2.001 

3.28 

0.75 

2.309 

3.49 

2.315 

3-65 

2.307 

3.64 

2.327 

3.68 

2.472 

4.07 

0.25 

1 .424 
2.08 

1 .417 
2.11 

1 .413 
2.10 

1  .425 
2.12 

1.484 

2.27 

1000 

9.0 

0.50 

4.510 

8.86 

4.4l6 

8.85 

*♦.397 

8.82 

4.452 

8.91 

4-73*4 

9-52 

0.75 

5.255 

10.88 

5-172 

10.86 

5.155 

10.83 

5.202 

10.91 

5-504 
11  .60 

0.25 

3.715 

6.79 

3.643 

6.70 

3-63*4 

6.69 

3-665 

6.73 

3.801 

7-01 

4000 

8.0 

0.50 

12.43 

25-12 

12.15 

24.82 

12.10 

24.74 

12.26 

24.96 

12.98 

26.31 

4000 

7.0 

0.50 

16.33 

31.4 

15.96 

31.1 

15.89 

31 .0 

16.13 

31.3 

17.21 

33-2 

7000 

**.5 

0.50 

54.01 

84.7 

52.53 

85-3 

52.15 

84.8 

53-39 

86.1 

58.48 

92.5 

0.75 

56.21 

97.2 

54.54 

98.2 

54.18 

97.8 

55-26 

98.8 

60.58 

105.1 

0.25 

48.83 

68.0 

W.52 

66 . 6 

*47.32 

66.4 

48.06 

67.1 

50.65 

70.4 

1 1 


Table  Caption 

Thermodynamic  results  for  binary  mixtures  of  exponential-six  fluids  representing  hydrogen- 
helium  mixtures  [11].  The  potential  parameters  are  1 1 1/kB=36.4K,  bn  =  3.43A,  an=11. 
1.  «12/kB=15.5K,  b12=3.37A,  «12=12.7.  «22/kB= 10.57K.  b22=2.97A.  <*22=13.6, 

where  the  subscript  1  denotes  hydrogen.  In  each  row  the  upper  results  are  pressures(GPa) 
and  the  lower  results  are  excess  internal  energies(kJ/mol).  Results  in  the  columns  headed 
Mixture  and  Ree  are  from  256-particle  Monte-Carlo  simulations  [11].  Results  in  the 
columns  headed  I.  II  and  III  were  obtained  by  using  Ross'  procedure  [12]  in  the 
perturbative  manner  described  in  the  text. 
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states  simulated  by  Ree  but  inferior  to  both  I  and  Ree's  empirical  rules  at  low 

temperatures .  We  do  not,  however,  regard  the  small  changes  in  accuracy  between  Ree,  I 

and  II  as  very  significant.  More  important  is  our  hope  that  our  approach  to  the  mixing 

rules  will  prove  useful  for  application  in  the  manner  indicated  above  to  non-conformai 

potentials  other  than  exponential-six . 
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g  «r)  =  5(a  r/b  )  (12) 

XX 

Substituting  these  relations  in  the  compressibility  and  energy  (or  virial)  equations  yields  the 
MRs 


III  <b3a  3>  ,  <«b3a  3e  a/(a-6)>  ,  <«b3a4/(a-6)> 

without  any  need  to  consider  only  the  part  of  the  potential  near  its  minimum. 
Unfortunately,  (11)  and  (12)  have  the  opposite  effect  to  what  is  expected  physically  at  high 
pressures:  they  indicate  a  principal  peak  in  g(ri  at  larger  r  for  softer  potentials  (i.e.  lower 
al. 

Ideally  one  would  like  to  test  the  MRs  proposed  above  by  simulations,  as  Ree  did, 
but  another  reliable  method  for  the  EPFs  is  Ross'  semi-empirical  variational  procedure 
[12],  which  we  have  used  to  obtain  the  tabulated  results .  Although  Ross '  procedure  is 
accurate  to  only  1-2%  in  absolute  terms  (and  seems  to  deteriorate  further  at  very  high 
pressures),  we  believe  that  it  can  accurately  distinguish  much  smaller  differences  between 
thermodynamic  properties  of  similar  potentials  8t  the  same  temperature  and  density.  Thus 
the  tabulated  results  are  based  on  the  assumption  that  the  change  in  exact  EPF 
thermodynamics  between  Ree's  MRs  and  I,  II  or  III  is  equal  to  the  corresponding  change  in 
the  Ross  thermodynamics . 

On  the  basis  of  these  results  we  conclude  that  III  gives  energies  and  pressures 
which  are  consistently  too  high,  I  is  (not  unexpectedly)  best  at  low  temperatures  but 
deteriorates  as  the  temperature  increases,  and  II  is  on  average  best  over  the  full  range  of 
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important  in  determining  the  energy  and  virial  equation  integrals.  We  therefore  use  in  the 
energy  (or  virial)  equation  the  Taylor  expansion  of  $  (r)  about  b  ^  and  demand  only  that  EPF 
and  mixture  results  obtained  by  integrating  the  first  two  non-vanishing  terms  of  the 
expansion  should  be  identical.  Apart  from  the  excluded  volume  MR  <b  >,  which  is  still 
obtained  exactly  as  before  from  consideration  of  the  compressibility  equations,  this 
procedure  leads  to  the  MRs 

I  <«b3>  ,  <«b3<x(a-7)/lct-6)> 
if  the  energy  equation  is  used  or 

II  <«b^a(o-7)/lot-6)>  ,  <«b^a(o^-56)/(o-6)> 

when  the  virial  equation  is  used. 

Rules  I  give  the  same  *  and  b  as  the  vdWl  rules  and,  since  the  a  are  all  in  the 
range  11-14,  an  a*  very  close  to  but  slightly  lower  than  the  value  obtained  by  Ree.  The 
functions  of  a  appearing  in  the  MRs  are  just  (♦lkl({=  1  ,a) |  ,  k=0,2,3,  where  <k)  denotes 
the  kth  derivative  with  respect  to  i .  A  straightforward  generalisation  when  there  is  more 
than  one  shape  parameter  would  be  to  obtain  the  extra  MRs  necessary  in  an  analogous  way 
from  higher  derivatives. 

For  exponential-six  potentials  it  might  seem  attractive,  from  a  purely  mathematical 
viewpoint,  to  assume  instead  of  (5)  and  (6)  the  scalings 

(1 11 


q  (r)  =  5(a  r/b  ) 
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3.  Generalisation 

One  generalisation  of  vdWI  theory  to  mixtures  of  conformal  diatomic  molecules  has 
recently  been  proposed  [8]  but  we  shall  be  concerned  here  only  with  non-conformai 
spherical  molecules.  There  are  at  least  two  reasons  why  one  may  wish  to  study  such 
mixtures.  First,  real  fluids  are  not  satisfactorily  represented  by  conformal  potentials  [9] 
and,  second,  the  effective  spherical  potentials  which  can  be  used  to  accurately  reproduce 
thermodynamics  of  some  non-spherical  molecules  are  not  in  general  conformal  even  when 
the  original  molecular  potentials  are  [10]. 

We  continue  to  scale  the  interactions  according  to  their  minima,  as  in  the  conformal 
case,  but  now  allow  for  differences  between  the  shapes  of  the  $  through  the 
introduction  of  a  set  of  parameters 

*  Ir)  =  <  *<  r/b  ,  a  )  (10) 

u  '1  *1 

The  simplest  and,  for  present  purposes,  the  most  relevant  realisation  of  (9)  is  the 

exponential-six  potential  where  there  is  only  a  single  parameter  o ^  for  each  ^  and  ♦<! , 

a)={6exp[a(1-{)]-ai_6}/(a-6).  Ree  [11]  has  recently  found  quite  empirically  that  for 

3  3  3 

exponentiai-six  mixtures  the  MRs  <b  >,  <«b  >  and  <«b  o>  lead  to  results  as  good  as 
those  of  vdWI  theory  for  Leonard- Jones  mixtures. 

We  wish  to  retain  the  assumptions  (5)  and  16),  which  are  physically  reasonable  at 
low  pressures,  but  observe  that  with  these  assumptions  and  an  exponential-six  pair 
potential  no  simple  scaling  of  the  energy  and  virial  integrals  occurs.  At  low  temperatures 
near  to  the  triple  point,  however,  we  expect  the  region  around  the  minimum  of  the 
potential  (coinciding  with  the  main  peak  of  the  radial  distribution  function)  to  be  most 
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Here  an  ambiguity  occurs:  if  A  is  considered  as  an  energy,  vdWl  yields  the  MR  <A>  but 
it  can  equally  well  be  considered  as  (potential  range)*’  leading  to  the  MR  <A3,V>.  The  latter 
rule  was  used  by  Evans  and  Hanley  [7]  for  v=12  and  gave  good  agreement  with 

simulations.  It  was  also  used  in  [4]  where  the  VdWl  results  for  inverse  power  potentials 
were  shown  to  be  quite  close  to  those  obtained  from  a  First  order  expansion  in  1/v  about 
the  Percus-Yevick  equation  of  state  far  binary  hard  sphere  mixtures. 
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term  to  g(r)  from  Eqs.  (5)  and  (6)  which  is  the  first  order  of  a  formally  exact  senes  in 
powers  of  fie: 

g(r)  =  Glr/I b)  {  1  *  fit  Hlr/b)  +  }  (7) 

where  H  is  assumed  to  be  'universal'  in  the  same  sense  as  G.  Using  (7),  Leland  et  a/. 

[2]  claim  to  obtain  both  of  Eqs.  (3)  and  (4)  from  consideration  of  the  energy  (or  virial) 

2  3 

equation  but  in  fact  it  seems  clear  that  one  obtains  Eq.  (4)  together  with  the  MR  <«  b  >. 

We  believe  that  a  more  satisfactory  way  of  obtaining  (3)  is  by  an  appeal  to 
thermodynamic  consistency.  If  we  substitute  (5)  and  (6)  in  the  compressibility  equation 

fi~'(3p/dp)p  =  1  +  p  l  x.x  J  dr  {  g(j(r)  -  1  }  (8) 

for  both  the  mixture  and  the  EPF  and  demand  that  both  give  the  same  result,  we  obtain 
Eq.  (3).  [We  remark,  incidentally,  that  substituting  (7)  in  the  compressibility  equations 
really  does  yield  both  (3)  and  (4)]  .  Thus  we  see  that  once  we  have  decided  to  represent  a 
conformal  mixture  by  a  conformal  EPF  and  made  the  assumptions  (5)  and  (6),  we  are 
naturally  driven  to  choose  the  vdWI  MRs  for  reasons  of  thermodynamic  consistency. 

For  pair  potentials  determined  by  a  single  parameter  only  one  MR  is  needed.  In  such 
cases  the  excluded  volume  MR  (3)  seems  to  be  more  important  than  the  energy  MR  (4). 
This  is  very  natural  for  hard  sphere  mixtures  where  the  potential  has  no  energy  scale  and 

3 

the  <d  >  MR,  which  for  hard  spheres  could  arise  from  consideration  of  either  the  virial  or 
compressibility  equations,  is  successful  [6].  Less  obviously,  it  also  seems  to  be  true  for 
inverse  power  law  potentials 

♦  (r)  =  A  lrv  =  <  (oJr)y  =  « Jo  Irf  ,  v>2  (9) 


2.  Justification 


Clearly  one  cannot  expect  to  'derive'  VdWl  theory  but  there  have  been  various 
attempts  to  understand  its  success.  Leland  et  al.  [M ]  note  that  it  arises  naturally  from 
the  vdW  equation  of  state  by  [hence  the  name  VdWl)  but  emphasize  that  its  applicability 
is  not  confined  to  fluids  satisfying  the  VdW  equation  of  state.  Smith  [5]  has  given  a  general 
perturbation  theory  based  on  the  reference  EPF  obtained  with  <embn>  and  <ePb^>  MRs. 
which  allows  systematic  corrections,  but  this  gives  no  insight  into  why  the  choice 
(m.n.p.q)-(0,3. 1 .3)  should  be  especially  accurate. 

Other  explanations  have  been  based  upon  the  assumption  that  the  suitably  scaled 
radial  distribution  functions  of  both  the  mixture  and  the  EPF  have  the  same  form,  i.e. 

g  (rl  =  G<r/b  I  (5) 

”  i  ij 


and 

g  (r)  =  Glr/b  ).  (6) 

where  the  function  G  can  depend  on  the  temperature  T=(kg/J)  1  and  the  total  number 
density  p  although  these  arguments  are  not  shown  explicitly.  If  relations  (1),  (2),  (5)  and 
(6)  are  substituted  in  the  energy  (or  virial)  equations  for  both  the  mixture  and  EPF,  then  the 
requirement  that  both  expressions  be  equal  leads  to  Eq.  (4)  [2,3].  This  seems  to  be  as 
satisfactory  a  justification  of  Eq.  (4)  as  one  can  hope  to  get,  but  the  arguments  used  by 
these  authors  for  Eq.  (3)  seem  flawed  by  comparison. 

Henderson  and  co-workers  [3]  merely  take  over  Eq.  (3)  from  the  similar  <d  >  MR 
Id  represents  sphere  diameter)  arising  from  consideration  of  the  virial  equation  for  hard 
sphere  mixtures  [6].  Leland  et  a/.  [2],  on  the  other  hand,  in  effect  use  a  correction 
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1.  Introduction 

Fluid  mixtures  are  of  great  practical  importance  but  theoretical  calculations  for  them 
are  generally  much  harder  than  for  single  spec  [1].  An  exception  is  mixtures  of 

conformal  spherical  molecules  not  too  different  in  size  having  pair  interactions  of  the  form 

4  Ir)  =  4  frr/b  )  (1) 

ij  ij  >j 

where  4>  is  a  universal  function  and  «  and  b  define  the  well  depth  and  the  separation  at 

*j 

its  minimum  respectively  for  the  interaction  between  species  i  and  j.  For  such  mixtures, 
the  van  der  Waals  one-fluid  (vdWI)  theory  [2]  is  believed  to  be  much  better  for 
predicting  thermodynamics  than  many  more  complicated  approximations  although  this  has 
admittedly  only  been  tested  for  Lennard- Jones  mixtures  with  ♦<{)={  -2{  [3] .  In 

vdWI  theory,  the  mixture  is  replaced  by  an  equivalent  pure  fluid  (EPF)  with  pair  interaction 

4  (r)  =  «  4Kr/b  ).  (2) 

X  XX 

Here  the  EPF  parameters  are  defined  by 

b  3  5  I  xxb  3  (3) 

X  '  I  'J 


and 


I 


I  x  x  *  b 

j  >]  ' 


(4) 


where  x  are  the  fractional  concentrations  of  the  components.  We  shall  have  occasion  to 

refer  to  many  different  mixing  rules  (MRs)  of  this  general  form  and  so  for  brevity  we 

3  3 

refer  to  13)  and  14)  in  an  obvious  notation  as  the  MRs  <b  >  and  <«b  >. 
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the  search  for  similar  procedures  for  non-spher ical  fluids  where  the 
situation  is  less  satisfactory  at  the  present  time,  although  a  certain 
amount  of  progress  has  been  made  in  recent  years  -  especially  for 
non-polar  one  component  diatormcs  [2],  It  is  the  purpose  of  this  note 
to  investigate  a  simple  scheme  for  mixtures  of  such  molecules. 


2.  Formulation  of  the  Problem 

The  interactions  in  a  n-component  mixture  of  symmetric  diatomic 
molecules  interacting  via  site-site  pair  potentials  have  the  form 

vij(R’B1*n2>  *ij(lR  *  1/2  V2]  *  1/2  V_22  l}  ’  ll) 

where  R  is  the  vector  between  the  centers  of  molecule  1  of  species  i 
and  molecule  2  of  species  j,  and  a  ,  k=l,2,  is  a  unit  vector  along 
the  axis  of  molecule  k.  £  .  and  Jt  are  the  bondlengths  of  the 
molecules  of  species  i  and  j  respectively,  and  ^(r)  is  the 
site-site  potential  between  a  site  of  species  i  ana  a  site  of  species 
j.  The  summation  in  Eq.  (1)  is  over  all  possible  +,  -  pairs  in  the 
argument  of 

We  shall  focus  here  on  the  particular  case 


e  i  J 
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-0-. 
-  y  j 


for  which  there  are  now  some  machine  computations  availuiie  for 
comparison  with  theory.  One  of  these  computations  ccrrespcrcs  to 
parameter  values  of  this  system  useo  for  the  moceii^c  or  a  lC~  - 

^2H6  mixture  around  zero  pressure  [3].  The  other  ccrrespcrcs  :c 

mixtures  of  molecules  having  the  same  c  arc  E  as  h, 

but  either  half  or  twice  its  bond  length  [5].  Tc  te  specif,c,  c.  :  *cr 
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N„  denoted  in  [51  by  N-N  is  0.32S2  while  t/o  =  G.164fc  anc  Ho  = 

£ 

C.6S84  respectively  for  the  molecules  denoted  in  [5]  by  N.h  anc  N — N. 

Ue  are  motivated  to  investigate  this  system  at  the  present  time 
by  a  desire  to  find  the  "best"  extension  of  the  recently  developed 
spherical i zation  procedure,  the  median  potential  [2, 6, 8, 9],  tc 
mixtures.  The  median  potential,  which  has  been  modified  in  various 
ways  [7],  has  been  found  to  give  good  results  for  single  component 
aiatomics  with  relatively  little  labor  [2,6-9].  The  analysis  is 
particularly  simple  and  accurate  for  one  component  L-J  molecules  where 
McGowan  [7]  found  that  the  best  modified  median  is  obtained  by 
splitting  the  potential  in  (2)  into  its  power  law  components.  This 
gives  a  median  with  a  universal  dependence  on  e,  a  and  l. 

There  are  various  ways  that  one  can  imagine  generalizing  the 
median  prescription  to  mixtures.  One  of  the  simplest  is  to 
spherical ize  all  angle  dependent  potentials  v.  about  the  centers 

J 

(or  some  other  suitable  points)  of  each  molecule  and  then  use  existing 

theories  of  spherical  fluid  mixtures.  This  is  the  path  followed  in 

[10]  for  "hard"  mixtures  of  dumbells,  spherocy 1 i nL ars  and  spheres. 

For  the  available  simulation  results  on  such  systems  the  resulting 

mixtures  of  "median"  hard  spheres  turned  out  to  be  very  nearly 

additive  -  for  which  the  F'ansoori-Carnahan-Starl  ing-lelanc 

approx  mat  ion  is  known  to  be  quite  accurate. 

Fcr  softer  interactions  there  is  no  simple  theory  even  for 

spherical  mixtures  except  for  the  case  of  van  cer  Waal s 1  type 

interactions  tf.-(r)  =  e  .  .  $(r/0..).  For  such  mixtures  the  "vcK 
TJ  ij 

1-fluid  thecry"  in  which  a  mixture  with  densities  ci  =  x.c 


represented  by  a  one  component  fluid  with  density  p  anc  potential 
e$(r/o) , 


-3  3 

o  =  >  x  .  x  .  o  .  .  , 

i  J  ij  ’ 

-  V*  3  .-3 

e  =  /  X  •  X  •  a..  e  .  .  /  a 

1  J  1 J  IJ 


(3) 

(4) 


is  at  least  as  accurate  as  much  more  elaborate  procedures.  This 
mixing  rule  has  been  extended  recently  by  Ree  [11]  to  mixtures  of 
particles  interacting  via  exponential-six  potentials  where  ^  .(r)  is 
of  the  form  e..<t>(r/c  ,  a  ),  4>(y,a)  =  (6  exp[a(l-y)]  - 

0  ■  o  *  J 

a/y  )/( a-6).  It  was  found  there,  by  trial  and  error,  that  adding  to 
(3)  and  (4)  the  rule 


a  -ZVj  eu“u  aij,(t  °3)  (5) 

gave  very  good  results  over  a  wide  range  of  densities  and 
temperatures.  This  is  encouraging  as  it  shows  again  that  suitable 
simple  schemes  do  work. 

Unfortunately  the  spherical  potentials  w.,(R)  obtained  b> 
"nedianizing"  the  angle  dependent  potentials  v.j(R,^,^)  are 
not  generally  given  by  a  simple  formula  for  which  (3)  and  (4)  can  be 
used  directly.  This  is  true  even  for  the  case  in  Eq.  (l)-(E).  It 
therefore  occurred  to  us  that  since  the  site-site  potentials  in  (2) 
are  of  simple  form  it  is  worth  trying  a  molecular  1-fluid  theory  in 
which  iT  and  T  would  still  be  given  by  (3)  and  (4),  and  It  would  be 
given  by  a  suitable  mixing  rule.  In  the  present  work  we  exclusively 
use  the  simple  rule 
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I  =£*,  i,  it) 

for  the  equivalent  bondlength. 

A  plausibility  argument  is  given  in  the  Appendix  for  a  more 
complicated  rule  (A. 16)  for  similar  to  (5),  but  for  all  mixtures 
where  simulations  are  available  the  £'s  obtained  from  (6)  and  (A. 16) 
are  almost  identical  (in  many  cases  exactly  so).  For  the  computer 
simulations  of  the  C0o  -  C0H.  mixture  treated  as  rigid  diatomics 

C  L  0 

[3],  the  two  lengths  £ ^  and  are  very  close  to  each  other  ana  so 
any  reasonable  mixing  rule  will  give  essentially  the  same  i.  For  the 
mixtures  of  [5],  on  the  other  hand,  one  is  just  testing  the  mixing 
rule  (6)  fori.  Unfortunately  the  number  of  available  results  is  not 
sufficient  to  permit,  at  this  time,  testing  of  different  possibilities 
for  T  and  indeed  there  appears  to  be  at  least  one  misprint  among  the 
published  data  in  [5].  Additional  computer  simulations  systematically 
varying  the  parameters  e--,  a..,  i,  and  x  would  be  most 

I  J  •  J  1 

useful . 


3 .  Results 

With  the  T,  1  and  1  given  by  Eqs.  (3),  (4),  and  (6),  w° 
sphericalize  using  MacGowan's  procedure  [7]  for  diatomic  lennard-Jcnes 
potentials  and  use  the  resulting  spherical  potential  in  Ross' 
semi-empirical  variational  procedure  [12]  to  obtain  the  Helmholtz  free 
energy,  and  from  it  by  numerical  differentiation  the  reported 
pressures.  We  also  calculate  the  excess  internal  energy  wherever  it 
is  available  from  the  simulation  .  We  remark  that  our  calculations 
are  very  quick,  requiring  less  than  one  second  of  CPU  time  on  a  DEC  20 
computer  per  thermodynamic  point. 
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Ue  compare  the  available  results  cf  molecular  dynamics 
simulations  of  CC-— C0h _[ 3]  with  our  calculations  in  Tables  I,  II 

C  C  c 

and  III.  Since  the  molecular  dynamics  calculations  for  CG^-C,^ 

have  as  one  of  their  purposes  to  obtain  the  set  of  o.  .,  e.  ., 

1 J  1 J 

i,j»l,2,  that  gives  the  best  intermolecular  potential  for  that  system; 
we  see  entries  in  Tables  I  through  III  corresponding  to  different 
values  of  those  parameters,  at  a  few  different  thermodynamic  points. 
Our  results  for  the  C0o-C„Hc  mixture  are  promising.  The 

C  C  0 

compressibility  factors  pV/MkT,  which  are  near  zero,  agree  with  the 
molecular  dynamics  calculations  within  about  C.2  which  correspcncs  tc 
differences  in  pressures  of  about  10  KPa,  and  the  excess  internal 
energies  per  molecule  in  units  of  kT,  U / NkT  to  within  0.2  which 
corresponds  to  differences  in  excess  internal  energy  of  about  C.5 
Kj/mol.  These  discrepancies  are  about  a  factor  of  2  larger  for  the 
pressures  and  about  the  same  for  U  than  those  observed  when  using  the 
median  procedure  for  the  pure  COj  and  fluids.  Also  the 

tendencies  for  our  calculated  values  are  the  same  as  for  the 
simulations. 

Ue  show  in  Table  IV  the  comparison  between  this  1-molecular 
fluid  approach  plus  "medi ani zat ion"  with  the  results  available  for 
pure  fluids  and  mixtures  of  molecules  of  elongations  i/c  of  0.6584, 
C.2292  and  0.1646,  respectively  [5].  As  alreaoy  mentioned  these 
molecules  are  denoted  in  [5]  by  N — N,  N-N  ano  IV. N  anc  all  have  the 
same  e  and  c  as  N„.  The  few  results  available  ao  not  allow  us  to 
draw  any  firm  conclusion,  he  believe  that  the  excess  internal  energy 
reported  in  [5]  for  the  pure  .\-N  fluid  is  in  error  (see  Table  IV). 
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Cur  results  show  that  the  approach  taken  coula  te  fruitful. 
However,  more  data  from  computer  simulations  are  needed  before  cne  can 
opt  on  a  firm  basis  for  one  of  the  possible  routes  to  an  extension  of 
the  vdW-l-fluid  approximation  now  in  use  for  mixtures  of  spherical 
molecules  to  mixtures  of  non-spherical  molecules. 
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Table  I 


urrber 

Ceos  i  ty 

p[ A~3 ] 

Ccncentra- 
ticn  of  CO^ 

X1 

t{  A) 

o22(A) 

o12(A) 

t12/k(°K) 

Ccirpress- 
i t i 1 i ty  Factor 
PV/NkT 

This  ND 

kork 

2.370* 

-0.54 

0 .01082 

0.50 

2.360* 

3.506 

3.248 

150.0 

-0.58 

-0.376  i 

2.349* 

-0.61 

C. 01082 

0.50 

2.360 

3.506 

3.300 

145.0 

-0.35 

C  .024 

C.C0865 

0.35 

2.356 

3.550 

3.300 

145.0 

-1.02 

-0.647 

C  .01037 

0.35 

2.356 

3.550 

3.300 

145.0 

-0.08 

C.074 

Comparison  between  the  molecular  dynamics  (MO)  calculations  [3] 

and  the  "med i an ized"  1-f  molecukr  approximation  explained  in  the  text 

for  a  mixture  of  CCL  (subindex  1)  -  C_H,  treated  as  rigid 

c  c  b 

o  ° 

Lennard-Jones  diatomics  with  =  2.370  A  and  i ^  =  2.349  A.  The 
compressibility  factors  are  near  zero  so  that  small  absolute 
differences  are  large  in  relative  terms.  The  temperature  was 

o 

r  =  241  X  on  =  2.989  A,  en/k  =  163.6  °K,  e  /k  =  137.5  °K. 


*Trese  three  values  are  < , ,  (>1+i„)/2  and  * „  respectively. 

lie  £ 

They  are  included  to  see  hew  cur  results  in  a  particular  case  are 
altered  tv  varying  7  frerr  t,  tc  Ln. 

1  c 
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Table  II 


Number  Concentra- 
Density  tion  of  CCU 

O  ^ 

r  y 


0.00833 
0.01302 
0.01023 
0 .01023 


e  12  /  k 


149.85 

149.85 

149.85 

120.0 


— 

Compress¬ 
ibility  Factor 
PV/NkT 

This  NO 

Work 

0.02 

0.C59 

-0.13 

0.066 

-0.54  - 

0.226 

-0.08 

0.099 

Excess  Internal 
per  Nolecule  in 
of  kT  L/NkT 
This 
Work 


Same  as  Table  I 
n..  =  3.3135  A. 


269.3  K,  =  3.035  A,  o  =  3.592  A, 
=  163.3  ° K .  e_/k  =  137.5  °K. 


Table  III 


C crr.pre ssibi  1  it>  Factor  Excess  Internal  Energy 

FV/NkT  per  Folecule  in  Units  of  kT 

U/NkT 

This  Work  N'D  This  Uork  NX 


0 .01369  1 

2.370 

0.33 

0.411 

-5.35 

-5 . 596 

C. 01303  1 

2.370 

-0.13 

C  .066 

- 

- 

0.01232  1 

2.370 

-0.48 

-0.186 

-4.79 

-5.028 

C. 01172  0.7 

2.364 

0.64 

0.723 

-4.53 

-4.818 

0.01021  0.7 

2.364 

-0.15 

0.129 

-3.88 

-4.166 

0.008696  C .  7 

2.364 

-u.42 

-C.102 

-3.21 

-3.599 

0.008914  '0.4229 

2.358 

-0.12 

0.124 

-3.60 

-3.698 

C. 008478  0.1656 

2.352 

-0.13 

0.127 

-3.68 

-4.090 

C. 009014  0 

2.349 

0.27 

0.445 

-4.63 

-4.800 

C. 008313  0 

2.349 

-0.20 

0.058 

- 

- 

C. 007812  0 

2.349 

-0.40 

-0.072 

-3.93 

-4.166 

C. 007608  0 

2.349 

-0.46 

-0.120 

-3.81 

-4.059 

Same  as  Table  I.  T 

=  269.2! 

3  X  The 

Lennard-Jones 

parameters 

are 

ojj  =  3.035  A,  =  3.520  A,  o^„  =  3.2775  A,  and  t^/k  = 

163.29  °K,  e22/k  =  137.49  X  e^/k  =  120  °K*  The  uncertaint>  1n 
the  compressibility  factors  is  estimated  to  te  +0.01. 
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Table  IV 


1  Temperature 

1  T 

1 

i 

System 

lie 

Compress i b i 1 i ty 
Factor 

PV/NkT 

This  Work  ND 

Excess  Internal  Energy 
per  Nclecule  in  Units 
of  kT  L/NkT 

This  Work  NT- 

102.6 

pure  N-N 

0.3292 

2.28 

1.91 

-5.89 

-3.67  + 

99.1 

pure  N — N 

0.6584 

19.47 

14.59 

-3.77 

-3.79 

1C1.C 

pure  N.N 

0.1646 

-1.66 

-0.92 

-6.60 

-6.22 

108.2 

equimolar 

N — N/N.N 

0.4150 

5.40 

4.45 

-5.13 

j 

-5.C3 

110.9 

equimolar 

N — N/ N-N 

0.4938 

9.05 

6.85 

-4 . 50 

-4.38 

'Comparison  between  NO  results  of  [5]  with  the  1-f  molecular  plus 


"medianization"  procedure  explained  in  the  text.  The  L-J  parameters 

o 

are  those  for  N~  [4],  a  =  3.31  A,  elk  =  37.3  K  for  all  mo.lecules. 

°_3 

The  density  for  all  cases  is  p  =  0.01779  A  .  As  explainec  in  text. 


the  it/ a  *  s  for  N.N,  N-N  ana  N — N  are  C.1646,  0.3292  and  C .  6584 


respectively . 


V<e  telieve  this  result  to  be  in  error. 


f\) 


REFERENCES 


I.  C.  Henaerson,  Annual  Review  cf  Physical  Chemistry,  25,  461 
(1974);  T.  W.  Leland,  P.  S.  Chappelear,  and  B.  U.  Gamson,  AIChE 
J.  8,  482  (1962);  T.  W.  Leland,  J.  S.  Rcwlinson,  and  G.  A. 

Sather,  Trans.  Faraday  Soc.  64,  1447  (1568).  For  a  recent 
review  on  fluid  mixtures  see  K.  E.  Gutbins,  K.  S.  Shing,  ana 
k.  B.  Street,  J.  Phys.  Chem.  87,  4573  (1983). 

C.  MacGowan,  E.  N.  kaisman,  J.  L.  Lebowitz  and  J.  K.  Percus,  J. 
Chem.  Phys.  80,  2719  (1584). 

3.  D.  Fincharr,  N .  Quirke  ana  D.  J.  Tilaesley,  to  appear.  ke  thank 
these  authors  for  making  the  results  of  their  calculation 
available  to  us  prior  to  publication. 

4.  P.S.Y.  Cheung  and  J.  G.  Powles,  Mol.  Phys.  30,  521  (1575). 

5.  K.  Nakanishi  and  H.  Tanaka,  Fluid  Phase  Equilibria  1_3,  371 
(1583). 

6.  J.  L.  Lebowitz  and  J.  K.  Percus,  J.  Chem.  Phys.  79,  443  (1983); 

N .  S.  Shaw,  J.  0.  Johnson,  and  B.  L.  Holian,  Phys.  Rev.  Lett. 

50,  1141  (1983) 

7.  0.  NacGowan,  to  be  publishea  in  J.  Chem.  Phys. 

8.  G.  0.  Willian.s,  J.  L.  Lebowitz  ana  J.  K.  Percus,  J.  Chen.  Phys. 

81,  2070  (1584). 

9.  C.  G.  Gray  and  C.  G.  Jcslin,  Chem.  Phys.  Lett.  101 ,  £48  (.  1 9  £  3 ) . 

10.  G.  0.  killiams,  J.  L.  Lebcwitz  ana  J.  K.  Percus,  to  te  publishea 

in  J.  Phys.  Chem. 

II.  F.  H.  Pee,  J.  Chem.  Phys.  _78,  409  (1983). 

12.  N.  Ross,  J.  Chem.  Phys.  71,  1567  (1979). 


38 


APPENDIX 


NIXING  RILES  FCR  HONGNUCLEAR  CONFORNAL  DIATCMCS 


For  ore-ccrriponent  hcmonuclear  diatomics,  the  energy, 
corrpress i b i  1  i ty  and  virial  equations  are  respectively 

Uex/N  =  2p J'dr  g000 ( r ) ^( r )  (A.l) 

kT(ac/ap)T  =  1  +  p  f dr  [gC0C(r)  -  1]  (A. 2) 

p  =  pkT  -  |  p  J'dr  [r  g0C0(r)  -  ( l/Sl)  g 1 CC (  r )  ]  aip/ar  (A. 3) 

with  the  following  analogous  results  for  [fixtures: 

Uex/N  =  2p 2xiXj  J*ar  g^° ( r )  i^(r)  (A. 4) 

kT(  3p  /  sp )  T  =  1  +  p  Exixj  /dr  Cg^j°  ( r )  -  1]  (A. 5) 


p  =  pkT  -  |  p  53xixj  Jar  [rg°j°(r)  -  ( X-i -//I)g -^0(r)]  d  ^,/cr 

(A. 6) 

here  p  is  the  total  number  density,  x.  the  fractional  concentra- 

r  rn 

tiers,  p..(r)  the  site-site  potential  energies,  a  ■  (r)  the 

1 J  ”  t  j 

1  GO 

site-site  racial  distribution  functions  with  g^,  the 
corresponding  1GG  spherical  harmonics  and 


where  i .  are  the  bcnolenaths. 
i 

ke  new  restrict  attention  to  site-site  potential  energies  having 
the  simple  scaling  property 

Vfr'  =  Ejj  -(rlc  ii)  'A.c! 

1  o  1  J  1  J 
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Clearly  it  is  possible  tc  obtain  a  sirrple  cne-fluic  rrcce  1  (analogous 
tc  van  cer  Waals  cne-fluic  for  atomic  licuics)  if  we  assume  that  the 
a.,  satisfy  the  simple  scalings 


000,  \  rCCC,  .  , 

g  (r)  =  G  (rio^) 


ICO/  X  r100 ,  ,  , 

g.  -  (r)  =  G  ( r/o  i  j  )  . 


(A.  1C 


loth  cf  these  may  be  included  in  the  single  statement 


j(r,n1.  1A)  =  G ( r / o i  _ ,  az) 


(A. 11 


for  the  full  angle-dependent  si te-centereo  correlation  function.  We 
assume  the  same  types  of  scaling  for  the  equivalent  one-component 
hcmcnuclear  diatomic  fluid,  that  is 


9  (  r ,  i'O  =  G  (r/o,  5^) 

Ur)  =7  'V(rlc) 


(A.  12 


(A. 13 


Inserting  (A. 12)  ana  (A. 12)  in  (A.l )  — (A.3)  and  (A. 8)  anr  (A.  11) 
ir  (A.4)-(A.6)  followed  by  eauating  (A.l)  with  (A. 4),  (A. 2)  with  { A .  £ 
anc  (A.3)  with  (A. 6)  yields  the  following  mixing  rules: 


where 


are  icertical  to  the  atomic  vck.  if  rules  arc. 


£ cs  .  r- .  '.u  ;  anc  ;  A . IE  ] 

■  r  fewest  i  rc  ly  ,  v,  A .  16 )  is  just  Ree's  rule  (5)  for  exponential-six 
ncratcmcs  with  a  replacea  by  x. 

*e  renark  that,  in  the  context  of  atonic  fluids,  (A. 14;  has 
usually  teen  justified  by  using  the  energy  equation,  or  al ternati vely 
the  virial  equation,  in  a  manner  analogous  to  that  described  above. 

On  the  ether  hanc,  the  argument  leading  frerr  the  compressibility 
equation  tc  (A. 15)  oces  not  seem  to  have  been  given  previously  for 
atomic  fluids.  Although  one  cannot  of  course  "derive"  approximate 
mixing  rules,  we  feel  that  our  use  of  the  compressibility  equation  is 
somewhat  more  satisfying  than  previous  approaches.  Once  we  decide  tc 
seek  an  equivalent  pure  fluid  and  make  the  scaling  assumption  for  the 
correlation  functions  then  reauiring  thermodynamic  consistency  leads 
naturally  tc  Eqs.  (A.  14)  and  (A. 15)  [and,  in  the  diatomic  case, 
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